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Abstract 

Multiparticle production in (2+1) dimensions is investigated. We show that 
in a small region around the threshold the perturbation theory becomes un- 
applicable due to infrared divergencies in a class of Feynman graphs with 
rescattering in final state. We develop a technique based on renormalization 
group for summing up leading logarithms and apply it to the (f)^ models with 
and without symmetry breaking and to the 0{N) cj)^ theory. 
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1. Multiparticle production in the conventional 0^ scalar theory is attracting considerable 
interest. Though the initial obserbvation [jl],0 was rather simple, only a few quantitative 
results has been obtained so far. Most of them concern the multiparticle threshold at which, 
for instance, formulas for the amplitude of production of n particle by an initial virtual one 
has been derived at the tree |]^-^ and one-loop levels |§-|§|. The tree amplitude exhibits 
a factorial dependence on the number of created bosons n. The correction from the first 
loop is of order An^ in comparison with the tree result (where A is the coupling constant) 
and one can expect the breakdown of the perturbation theory at large n. There have been 
several attempts to take into account contributions from all orders of the perturbation theory 



T0[| . In particular, the amplitude has been calculated at leading order in in the 0{N) 



symmetric scalar theory at large N [|T0[. However, it remains unclear whether the 1/N 
expansion is applicable at large number of final particles n. 

In this paper we consider the multiparticle production in (2+1) dimensions, where beside 
the above-mentioned inapplicability of the perturbation theory at large n there is another 
problem related with the infrared divergence of a certain kind of loop graphs which breaks the 
conventional pertubative expansion in a region close enough to the threshold. The graphs 
of this kind contains loops related to the rescattering among final soft particles. It is a 
pecular feature of (2+1) dimensions that right at the threshold, rescattering graphs diverge 
logarithmically, so even at small number of final particles n calculation of the amplitude 
around the threshold requires a nontrivial summing an infinite set of graphs. 

We find a suitable technique to do such kind of summing which is a modification of 
the conventional renormalization group (RG). This technique allow us to sum up leading 
logarithms from all orders of the perturbation theory. By this technique we obtain the 
amplitude 1 ^ n in the (f)'^ models with both broken and unbroken reflection symmetry, and 
in the 0{N) 0^ theory. We check the leading in 1/A^ amplitude in the latter model and 
show explicitly that the 1/N expansion breaks down at large number of final particles. 

2. Let us consider the 0"^ theory in (2+1) dimensions. 
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L = \{d,ct.f-\ct^'~\ct.\ 

where, for convenience, we suppose that the mass of the particle is equal to unity. This is 
the theory of interacting relativistic bosons. 

To describe bosons in the low-energy limit, one writes the following effective lagrangian 
in term of a non-relativistic bosonic field \E', 

where g is some, yet to be determined, effective coupling 

Non-relativistic bosons, interacting via a delta-like potential (as in eq. (P), have been 
known for a long time as an example of a non-relativistic system with dimensional transmu- 



tation ||Tl[. In fact, from dimension counting in (2+l)d (in non-relativistic kinematics the 
relation between the dimensions of energy and momentum is [E] = [p]'^) it is implied that 
vl>t\i>tvi>vl> is an marginal operator and the coupling constant g is dimensionless. Appearently, 
there is no scale parameter in the theory described by the lagrangian (|l]). However, this 
scale exists and is simply the boson mass (which is unity in our notation), which plays the 
role of the "ultraviolet" cutoff in the effective theory. 

To make contact between the effective lagrangian and the initial Lorentz-invariant one, 
one compares formulas for the amplitude of elastic scattering of two bosons computed in 
both theories. This results the following relation between g and A, 

9 = 3A/8. (2) 

Considering g as the bare coupling, or the coupling at the scale of ultraviolet cutoff, one can 
question about the evolution of g as a function of the momentum scale. For this purpose we 
introduce the running coupling constant g{t), which has the physical meaning of the strength 
of the interaction between bosons at the momentum scale p = e~*, {t = ln{l/p)). There is 
only one Feymnan diagram (shown in fig. |l]) that makes contribution to the corresponding 
beta function. Simple calculations yield the following RG equation. 



dg{t) 



dt 



TT 



which has the solution 



9{t) 



1 + H 

TT 
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) 



(3) 



where we have made use of the initial condition on g{t), g{^) = g, and g is defined by eq. 
(0). From eq. (^) one sees that the intensity of the interaction between bosons decreases as 
the momenta of particles tend to zero. Later, we will demonstrate that this property holds 
for a more general case of 0{N) cf)^ model, but in the theory with symmetry breaking the 
behavior of g{t) is diametrically opposite. 

The coupling constant changes considerably form its initial value only when the momen- 
tum scale is exponentially small on A, t = ln(l/p) ~ 1/A. So, the renormalization group is 
suitable for considering the regime A 0, At ~ 1. All further considerations will be done 
for this particular regime. 

In fact the fiow of the effective coupling, eq. can be obtained by a more trivial method 
by direct summing bubble graphs: in non-relativistic theories only these graphs contribute 
to the elastic scattering 2-^2. However, in calculation of the amplitude of multiparticle 
production around the threshold more complicated diagrams are involved and the problem 
cannot be reduced to the summing of bubble graphs. For example, for the 1^3 process the 
one- and two-loop diagrams that make contribution in our interesting regime is presented 
in fig. 0. To deal with these processes the RG technique is required. On the diagramatic 
language, the renormalization group corresponds to summing leading logarithms, i.e. terms 
proportional to (At)" from the whole perturbative series. 

Let us try to describe the production of n final particles in term of the non-relativistic 
bosonic creation operator \E'''". The only relevant candidate is An^'^^. The following relation 
should take place. 



(n|</)|0) = (n|A„^tn|o) 



eff, 
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where the matrix element in the l.h.s. is written in the initial 0^ model and the r.h.s. is 
understood as a matrix element in the effective non-relativistic theory. Form this equation 
one finds that An is equal to the 1 ^ n amplitude when the spatial momenta of final particles 
is small, but not exponentially small (so the logarithms do not appear in loops). Let us for 
simplicity consider small enough n where An coincides with the tree 1 ^ n amplitude , 

/^ \ ("-i)/2 

An = A'r = n\ i^-j . (4) 

When characteristic scale of momenta of final particles is exponential, the r.h.s. of eq. 
(1) is substantially renormalized by loops. One can treat this renormalization by introducing 
the running vertex An{t) and solving the corresponding RG equation, 

= ^ g{t)An{t). (5) 

at Z TT 

Graphically, this equation is represented in fig. ^. Substituing the function g(t) found in 
to eq. (H), the function A(t) can be found easily, 

An{t)=An[l + '-ftj ' =^T{l + ^^tj ' . (6) 

So, we have found that the RG techniqe allows us to calculate the 1 — n amplitude An{t) 
in a region close to the threshold where At ~ 1. Note that exactly at the threshold, i.e at 
t = +00, the amplitude vanishes. 

3. All the said above is equally applicable for the case of broken symmetry, 

with the only exception that for calculating the bare coupling g entering into the lagrangian 
of the effective theory from the initial model one must take into account not only the 
diagram with a four-boson vertex, but also diagrams with two three-boson vertices, i.e. 
those with exchange of a virtual particle in s-, t- or u-channels. In contrast with the case 
without symmetry breaking, the resulting amplitude is negative, which means the attractive 
character of the force between boson at low energies. One has 



3A 

(we assume v"^ = (2A)^^, so the mass of the boson is unity). Formula (Q) imphes that the 
strength of the interaction increases with t, 



9it) = T^ = --[l--t] . (7) 



Litterally, eq.(|^) predict an infinite couphng constant at an exponentially small momentum 
scale of 

Po = exp(— 27r/3A). 

This fact is a direct analog of the Landau pole in field theories without asymtotic freedom. In 
reality, it is a manifestation of the existence of a two-particle bound state in our model (recall 
that at least one bound state exists in every, arbitrarily weak, two-dimensional attractive 
potential). One can show that in our case, the energy of the bound state is of order pg. 
The 1 — > n amplitude is now 

n(n — l) 

27r 



Mt) = AT fi - ■ (8) 



It is interesting to compare our result for the case of large n, n ^ 1 to that obtained in ref. 
0] for (2+1) dimension by a different method. In the case when At is small. At -C 1, our 
formula in fact reproduces the result of 

However, if At is comparable with 1, the formula (P) does not coincide with that obtained 
in ref. . We consider it as a counter-argument to the claim of ref. [§] . 

4. Application of the technique described above the multi-component model requires 
a slight modification. The lagrangian of the 0{N) cf)^ model, 

1 1 A 

L = ^{d^(l)a)idM - ■:^<Pa<Pa " ^(</'a</'a)^ (9) 



6 



(a = 1 ... is the isospin index) contains one coupling constant, A. However, if one tries to 
write down the most general 0{N) symmetric effective non-relativistic lagrangian, one sees 
that there may be two different effective couplings, gi and g2, which correspond to the two 
possible structures of the potential term, 

gi is the low-energy elastic scattering scattering amplitude, singlet in the s-channel, while 
g2 determines amplitudes, singlet in t- and u-channels. From the initial lagrangian (|^) one 
obtains the bare value of the coupling constants, 

91 = 92 = ^- (11) 

The fact that the bare values of gi and g2 are equal is the remnant of the crossing symmetry 
of our intial lagrangian (^). However, the evolution equations for the running coupling 
constants are rather different, 

^ = -[N9!it)+Ag,it)g2m (12) 
at Tx 

= ^-nUt). (13) 

at 71 

In order to simplify the RG equations, let us introduce instead of gi{t) a linear combination 
of the two couplings, 

9it) = giit) + ^g^it). 
The RG equation for g(t) is simplier than that of gi{t), eq. ([12|) , 

MO = Hg^t). (14) 

at 71 

The solution to the RG equations, eqs. ([I3| ) and ([l^ , which satisfies the initial condition 
(0), can be found. 



In analogy with the simple (p^ case, the production of n soft bosons from an initial particle 
with isospin a can be described by an effective operator y4„,\E']j(\l'||\E'||)*-"~^''/^ {n must be odd). 
After some calculations we obtain the following RG equation for An{t), 



dt TT 



^(iV + n - l)g(t) + (n - 1)2 (l - 1 ) g,(t) 



An{t). 



(15) 



Having substitued the formulas for (72 (^) and g{t) to eq. ([TsD and solved it, one obtains the 
dependence of the 1 — n amplitude on the logarithm of the characteristic momentum of 
final particles, ln(l/p) = t 



-^{N+n~l) 



Air 



2N 



(N-l) 



(16) 



Recall that the regime where we are working is A 0, At ~ 1, provided other parameters 
as n and are fixed. 



Let us consider the large limit and show that the result of ref. can be reobtained 
in (2+1) dimensions by our technique. Consider the case when n much smaller than both 



the number of boson flavors and the inverse coupling constant, eq.(|l^) reduces to a much 
simplier formula. 



^ NX 

1 + 1 

Stt 



-(n-l)/2 



(17) 



Remember that the coupling constant enters into the tree amplitude by the factor of A^" ^•'Z^, 
one can rewrite eq.(0) into the following form. 



(18) 



where 



A; 



A 



is just the renormalized singlet scattering amplitude. One sees that for a small number of 
final particles, n <^ N, leading order in 1/A^ result is given just by the tree-level formula 
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where the couphng constant is replaced by the renormahzed one. This is precisely the result 



of [0 in the particular case of (2+1) dimensions. 

However, formula (0) has been deduced provided that n N. If the number of final 
particles is comparable with A^, the effect of loops is obviously not a simple renormalization 
of the coupling constant. One can find from (|TBp that the correction to the leading order 
(|17|) is proportional to /N . When the number of final particles becomes comparable to 
the number of their sorts, the expansion becomes totally unreliable. One can expect 
that the breakdown of the expansion is not a pecular feature of (2+1) dimensions but 
holds also in the realistic (3+1) dimensional theory. 

6. So, we see that the renormalization group is a poweful mean for investigating the mul- 
tiparticles amplitudes in (2+1) dimensional scalar field theory at and around the threshold. 
The exact formula for the amplitude, if ever be found, must incorporate the information 
obtained here from the renormalization group equations. 
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FIGURES 

FIG. 1. The only graph contributing to the lowest-order beta function in the effective 

non-relativistic theory. 

FIG. 2. One- and two-loops Feynman graphs that contribute to the 1^3 amplitude in the 
interesting regime. 

FIG. 3. A graph that make contribution to RG equation for 1 — > n amplitude. There are 
n{n — l)/2 diagrams of this type, corresponding to different possiblities to choose of a pair of final 
particle to rescatter. 
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